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Abstract 

We study the drift of slow variables in a slow-fast Hamiltonian system 
with several fast and slow degrees of freedom. For any periodic trajectory 
of the fast subsystem with the frozen slow variables we define an action. 
For a family of periodic orbits, the action is a scalar function of the slow 
variables and can be considered as a Hamiltonian function which generates 
some slow dynamics. These dynamics depend on the family of periodic 
orbits. 

Assuming the fast system with the frozen slow variables has a pair of 
hyperbolic periodic orbits connected by two transversal heteroclinic tra- 
jectories, we prove that for any path composed of a finite sequence of slow 
trajectories generated by action Hamiltonians, there is a trajectory of the 
full system whose slow component shadows the path. 

1 Introduction 

We consider a slow-fast Hamiltonian system described by a smooth Hamiltonian 
function H (p, q, v, u; e) . This system is slow-fast due to a small parameter in the 
symplectic form 

f2 = dp A dq H — dv A du . 

e 
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Therefore the equations of motion take the form 

dH dH 



dp dq m 

dH dH y } 

u = e— ) v = -e— . 
dv du 

Equations of this form often arise after rescaling a part of the variables in a 
Hamiltonian system with the standard symplectic form. 

The variable (p, q) are called fast and (v, u) are slow. We assume that the 
system has m + d degrees of freedom, where m is the number of fast degrees of 
freedom and d is the number of slow ones. 

After substituting e = into equation ([!]) we see that the values of (v, u) 
remain constant in time and the system can be interpreted as a family of Hamil- 
tonian systems with m degrees of freedom which depends on 2d parameters. We 
call it a frozen system: 

. _dH dH 

^ dp ^ dq (2) 

ii = 0, i) = 0. 

The case when the fast system has one degree of freedom is relatively well 
understood. Indeed, in this case the frozen system typically represents a fast 
oscillator. The averaging method can be used to eliminate the dependence on the 
fast oscillations from the slow system. Therefore trajectories of the slow system 
are close to trajectories of an autonomous system with 2d degrees of freedom over 
very long time intervals (see e.g. [1, 2j). 

Points of equilibria of the frozen system form surfaces called slow manifolds. 
Normally hyperbolic slow manifolds persists and normally elliptic slow manifolds 
do not in general. In both cases the dynamics in a neighbourhood of a slow 
manifold can be described using normal forms (for a discussion see e.g. [5]). 

The case when the fast system has more than one degree of freedom is notably 
more difficult. The effect of the fast system on the slow variables strongly depends 
on the dynamics of the fast system. If the frozen fast system oscillates with a 
constant vector of frequencies, generalisations of the averaging method can be 
used [8j [9]. The averaging method can be also used if the frozen system is 
uniformly hyperbolic [I] or, more generally, if the frozen system is ergodic and 
time averages converge sufficiently fast to space averages (?]. In all these cases 
the dynamics of the slow variables is described, in the leading order, by the vector 
field obtained by taking an average of the slow component of ([I]) over the space 
of fast variables 

/dH\ /dH\ 

This approximation strongly relies on the fact that in an ergodic system the time 
average over a trajectory equals the space average for almost all trajectories. 
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The approximation error strongly depends on the rate of convergence for time 
averages. If the number of fast degrees of freedom is larger than one, there is 
no reason to expect that the time average over a periodic orbit converges to the 
average over the space. Therefore we should expect that the slow component of a 
trajectory whose fast component stays near a periodic orbit of the frozen system 
may strongly deviate from the averaged dynamics described by ([3]). Moreover, 
we note that periodic orbits are dense in the case of an Anosov system. 

In this paper we assume that the frozen system has a compact invariant set 
bearing chaotic dynamics of horseshoe type created by transversal heteroclinics 
between two saddle periodic orbits. This situation typically arises when a periodic 
orbit has a transversal homoclinic. 

In this invariant set hyperbolic periodic orbits are dense and every two periodic 
orbits are connected by a heteroclinic orbit. We select a finite subset of periodic 
orbits with relatively short periods. We construct trajectories of the full system 
which switch between neighbourhoods of the periodic orbits in a prescribed way. 
We show that the slow component of such trajectories drifts in a way quite similar 
to trajectories of a random Hamiltonian dynamical system with 2d degrees of 
freedom. 

The trajectories constructed in this paper strongly deviate from the averaged 
dynamics. We think this mechanism is responsible for the largest possible rates 
of deviation. 

A similar construction is used in [6] for studying drift of the enrgy in a Hamil- 
tonian system which depends on time explicitly and slowly. In particular, it 
was shown in [6] that switching between fast periodic orbits indeed provides the 
fastest rate of energy growth in several situations. 

The rest of the paper has the following structure. In Section [2] we state our 
main theorem and discuss its application to systems with one slow degree of 
freedom. In Section [3] we describe slow dynamics of the full system ([1]) near 
a family of periodic orbits of the frozen system. The description is based on 
an action associated with the frozen periodic orbits and can be of independent 
interest. The central ingredient of the proof of the main theorem is preservation of 
normally hyperbolic manifolds formed by families of uniformly hyperbolic orbits 
of the frozen system which is explained in Section HI In this section we explain 
how symbolic dynamics can be used to describe the dynamics of the full system 
restricted to an invariant subset close to the hyperbolic invariant set of the frozen 
system. The discussion is based on ideas of jS] • Section 5 analyses the long time 
behaviour of the slow component of the full dynamics. The last section of the 
paper finishes the proof of the main theorem. 
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2 Accessibility and drift of slow variables 



The total energy is preserved, so we study the dynamics on a single energy level. 
Without any loss in generality we may consider the dynamics in the zero energy 
level 

M £ = {H(p,q,v,u;e) = 0} . 

First we state our assumptions on the dynamics of the frozen system. Let D C M. 2d 
be in a bounded domain. We assume 

[Al] the frozen system has two smooth families of hyperbolic periodic orbits 
L c {y, u) C jMo defined for all (v, u) G D, c G {a,b}. 

[A2] the frozen system has two smooth families of transversal heteroclinic orbits: 

r ab (v,u) C W u (L a (v,u))nW s (L b (v,u)), 

T ba (v,u) C W u (L b (v,u))C)W s (L a (v,u)), V(t> , u) £ D . 

We note that under these assumptions the frozen system has a family of uniformly 
hyperbolic invariant transitive sets Ar VjV ,\, also known as Smale horseshoes. For 
every (y, u) G D, this set contains a countable number of saddle periodic orbits, 
which are dense in Ar VjU y Moreover, every two periodic orbits in A(„ )lt ) are con- 
nected by a transversal heteroclinic orbit, which also belongs to At V)U y It is well 
known that the dynamics on the Smale horseshoe can be described using the 
language of Symbolic Dynamics. We define 

A:= |J A {v>u) . 
(v,u)eD 

Before stating our main theorem we give a couple of definitions. 
Definition 1 The action of a periodic orbit L c is defined by the integral 



J c (v,u) := (b pdq. 

J L c (v,u) 

The function J c is independent of the fast variables and can be considered itself as 
a Hamiltonian function which generates some slow dynamics of (v, u) variables: 

/ 1 dJ c , 1 dJ c 

V= -T C &I> M= 7^7' (4) 

where ' stands for the derivative with respect to the slow time r = et, and 
T c is the period of L c . System (jlj) is Hamiltonian with the symplectic form 
oj c = T c (v, u)dv A du. Alternatively the equations can be interpreted as a result 
of a time scaling in a standard Hamiltonian system. 
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In the next sections we will show that for properly chosen initial conditions the 
slow component of the corresponding trajectory of ([1]) oscillates near a trajectory 
of this slow Hamiltonian system. 

Inside the Smale horseshoe there are infinitely many periodic orbits connected 
by transversal heteroclinics. Each periodic orbit has an action associated with 
it. We select a finite subset of periodic orbits and consider the collection of their 
actions. In general we should expect all those actions to be different. 

In this paper we prove that there are trajectories of the full system such 
that their slow components follow any finite path composed of segments of slow 
trajectories generated by actions. Those trajectories of the full system shadow 
a chain composed of the periodic orbits and heteroclinic trajectories and spend 
most of the time near periodic orbits of the frozen system. 

Let us give a definition of an accessible path and then state the theorem. 
Consider a finite family of functions J k '■ D — > R, k = l,...,n. Let be 
the Hamiltonian flow with Hamiltonian function J k and the symplectic form 

— Tk(v, u)dv A du where T k > is the period of the corresponding orbit. For 
every point z = (v, u) G D we define 

a k (z) = sup{ r : &k(z) G D for all r' G (0, r) } , 

which is the time required to leave the domain D. If the trajectory is defined 
for all r > we set cr k (z) = +oo. Obviously, a k (z) > for any z and k due to 
openness of D. 

Definition 2 We say that T : [0,T] — > D is an accessible path if T is a piece- 
wise smooth curve composed from a finite number of forward trajectories of the 
Hamiltonian systems generated by J k . 

More formally, T is an accessible path if there are = To < T\ < . . . < rjv = T 
such that the sequence of points := r(rj) breaks the curve T into trajectories, 
i.e., for every i < N, there is ki, 1 < ki < n, such that for r G fa, r i+ i] 

r(r) 

Of course, the curve T is well defined only if 

< T i+ i ~ Ti < CTkiiZi) 

which ensures that the trajectories do not leave the domain D. 

Theorem 1 If D is a bounded domain in M? d , the frozen fast system satisfies 
assumptions [Al] and [A2], {Jk} k=i is a set of actions corresponding to a finite 
set of frozen periodic orbits in A, and T is an accessible path, then there is a 
constant C > and e > such that for every e < e there is a trajectory of the 
full system (TjQ) such that its slow component z(t) satisfies 

||z(t)-r(et)|| <C e 

provided < t < E~ l T . 
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Definition 3 For any zo,z\ G D, we say that z\ is accessible from zq via the 
system {Jk} if there is an accessible path such that T(0) = Zq and T(T) = z\. 

In the case d — 1 the accessibility property has a simple geometrical meaning 
since trajectories of the Hamiltonian systems generated by Jk are level lines of the 
functions Jk- In this case the theorem provides trajectories which follow segments 
of the level lines. The main obstacle for the drift in the slow space is provided 
by level lines common for all Jj.. 

Corollary 1 Consider actions generated by two periodic orbits, a and b. Those 
level lines of J a ^, which are inside D, are closed curves. The non-singular level 
lines form rings (or disks), D a and D^. Let V = D a fl A> C D. If J a and Jf, do 
not have common level lines, then any point z\ G V is accessible from any point 
z G V. 

Corollary 2 Under the same assumptions. Let us take any finite family of open 
sets Vi C V , which do not depend on e. Then for all sufficiently small e, there is 
a trajectory which visits all the sets Vi. 

If the energy set M. £ is compact the slow dynamics never leaves a bounded 
set. If at the same time D is a connected set, natural questions arise: Is there 
a point in D which is not accessible from every other point in D? Is there a 
trajectory such that its slow component is dense in D7 

3 Actions and first return maps near periodic 
orbits of the frozen system 

Now consider the cylinder formed by periodic orbits of the frozen system: 

S c ,o= |J L c (v,u)cM . (5) 

(v,u)eD 

Let 7 £ denote a trajectory of the full system jTJ and n s : ]R 2m+2d — > R 2d the 
projection on the slow variables. 

In the next section we will prove that some trajectories stay in a neighbour- 
hood of S Ct o for a very long time and provide a detailed description for them. In 
this section we show that in this case the evolution of the slow component tt s ^ £ 
approximately follows a trajectory of the slow Hamiltonian flow generated by 
the action J c . 

Lemma 1 Let L c be a family of periodic orbits of the frozen system. If j £ is a 
family of solutions of the full system (JT]) such that 
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(i) there are z$ G D and Cq > such that 

||7r s 7 e (0) - zq\\ < C e, (6) 

(ii) there are constants C\ > and tq < a c (zo) such that 

dist( 7£ (t),5 c , ) < Cte Vt G [0,6-Vo], (7) 

i/ien t/iere is C 2 > suc/i that 

||7r. 7 .(f) -*?(*>) II <C 2 s (8) 

/or all t e [0,e' l T ]. 

Proof. We write (p c (t, v, u), q c (t, v, u)) to denote a periodic solution of the frozen 
system and use T c (v, u) for the corresponding period: 

p c (t + T c (v,u),v,u) = p c (t,v,u), , g . 
q c (t + T c (v,u),v,u) = q c (t,v,u). 

Then the action of the periodic orbit L c is given by the following integral 

J c (v,u)= [ TC p c ^dt. (10) 



dt 

Since L c belongs to the zero energy level we have a useful identity: 

H(p c (t,v,u),q c (t,v,u),v,u;0) = (11) 

for all (v,u) e D and all t e R. 

Let S denote a smooth hypersurface in K 2m+M transversal to the flow of the 
frozen system such that every periodic orbit of the family L c has exactly one 
intersection with E. Let Mj = 7 e (tj) be a sequence of consecutive intersections 
of 7 e with E and consider the slow components of those points: Zi := 7r s Mj. 

We note that inequality © and the smooth dependence of p c (s, z), q c (s, z) on 
z imply that there is C3 > such that for every % there is s% such that 

\\Mi - (p c (si,Zi),q c (si,Zi),Zi)\\ < C x e. 

Since solutions of differential equations depend smoothly on the initial conditions 
and vector field, the segment of %(t), U < t < t i+ i is close to L c {z^)\ 

%(t) = (p c (si + t-U, zi), q c (si + t-U, Zi),Zi) + 0(e) (12) 

and the time of the first return to the section E is close to the period of the frozen 
trajectory: 

t, i+l -t i = T c {z i ) + 0{e). (13) 
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Now we estimate the displacement Zi+i — z,{. We write Zi = (v, u) and z i+ i = (v, u) 
to shorten the notation. Integrating the slow component of the vector field along 
the exact trajectory and using (DD) we conclude 



u — u 



v — V 



ti+1 



udt 



ti+l 



vdt 



Tc(v,u) QJJ 

p c (t,v,u),q c (t,v,u),v,u 
To(v,u) Q£[ 



dt + 0(e 2 



du 



(14) 



dt + 0(e 2 



p c (t,v,u),q c (t,v,u),v,u 



where the error terms come from replacing the exact trajectory by the frozen one 
and from the difference in the return time, see (TT2|) and (TXHT) . The integrals in the 
right hand side can be expressed in terms of derivatives of the action defined by 
integral ({TO]) . Indeed, differentiating (fTOj) with respect to u, integrating by parts 
and taking into account we get 



dJc 
du 



dp c dq c dq c dp, 



n du dt du dt 
1 'dpcdH dq c dH 



o 



du dp du dq 
Then differentiating identity ( fill we get 



* 1 dt 

dt 



dp c dH dq c dH 



dH 

du 



du dp du dq 

where the derivatives are evaluated at (p c (t, v, u), q c (t, v, u), v, u). Consequently 



dJc 
du 



du 



dt 



p c (t,v,u),q c (t,v,u),v,u 

Repeating these arguments with u replaced by v we also get 



dJc 
dv 



T C (V,U) QJJ 

dv 



dt. 



p c (t,v,u),q c (t,v,u),v,u 

Substituting the last two equalities into (TBI) we arrive to 



dJ 

u = u-e—^ + 0(e 2 ), 
dv 



dJ 

v = v + e— - + 0(e 2 ) 
du 



(15) 



We see that the displacement between two consecutive intersections of 7 e with 
section E is approximated by the time-eT c shift along a trajectory of the Hamil- 
tonian vector field (j3J) generated by the Hamiltonian function J c : 



i m = ^ + $f c (5 4 ) + 0(£ 2 ). 
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Inequality © implies that z = z + 0(e). Then a rather standard stability 
estimate can be used to show 

z { = & c £Tc (z ) + 0(e) for < i < e'W (16) 

Finally, we note that 2i = TT s le(ti) where = iT c + 0(ei) due to (fl3l) . Be- 
tween intersections with E the slow component 7r s 7 e changes by a value of the 
order of 0(e). Therefore the estimate is extendable to values of t between the 
intersections and (JHJ) follows immediately. □ 

We note that J c is preserved by and therefore J c o n s is an adiabatic 
invariant for the restriction of the full dynamics on a neighbourhood of S Ct Q. 

In general, we do not expect the estimates to be valid on time intervals longer 
than stated by Lemma [U For example, note that a trajectory of J c may leave 
the domain D in finite time. 

It is interesting that under additional assumptions J c (t) := Jc^sle^)) may 
stay near its initial value, J c (0), for much longer time. 

Indeed, consider the case of one slow degree of freedom, d — 1, and suppose 
that level lines of J c are closed curves on the plane of (v,u) variables. Then 
<y c (z) = +oo for all z G D. In the next section we will show that the full system 
has a locally invariant cylinder S CjE close to S C £. Then equations CH]) describe 
a Poincare map on the section defined by intersection of S cx and S. In the 
case of one slow degree of freedom we may suppose that the map (THI) satisfies 
assumptions of the KAM theorem, then J c (t) will stay close to its initial value 
forever. Indeed, under these assumptions the trajectories on S Ct£ are trapped 
between two KAM tori. 

We also note that averaging theory can be used to study the dynamics on S C)£ 
for d > 1. 



4 Dynamics of the frozen system and normal 
hyperbolicity 

The following arguments are based on [6]. Let w = (p,q) and z = (v,u) to 
shorten notation. Then the frozen system (TSJ) has the form 

w = G(w,z), (17) 

where G is expressed in terms of partial derivatives of H for e = 0. The Hamil- 
tonian function H(w, z) is an integral of system (fl71) . 

Let system fflTj) have two smooth families of saddle periodic orbits L a : w = 
w a (t, z) and Lf, : w = Wb(t, z) for all z G D. Assume that both families belong to 
the zero energy level J^Aq. Take a pair of smooth cross-sections, S a and S^, which 
are transverse to the vector field and such that each periodic trajectory L a (z) 
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and L b (z) has exactly one point of intersection with the corresponding section. 
Denote the Poincare map on S c near L c as U cc (c = a, b). The Poincare map is 
smooth and depends smoothly on z. 

We assume that for all z £ D the frozen system has a pair of transversal 
heteroclinic orbits: r afe C W u (L a ) n W s (L b ) and r fea C W u (L b ) n W s (L a ). 

Let U ab and be maps defined on subsets of S a and S b by following orbits 
close to T ab and T ba , respectively. Therefore U ab acts from some open set in S a 
into an open set in while acts from an open set in E& into an open set in 
S a . There is a certain freedom in the definition of the maps U ab and ITb a . Indeed, 
each of these maps acts from a neighbourhood of one point of a heteroclinic orbit 
to a neighbourhood of another point of the same orbit, therefore different choices 
of the points lead to different maps. 




Figure 1: Poincare maps near two periodic orbits 

When the maps are fixed, every orbit that lies entirely in a sufficiently small 
neighbourhood of the heteroclinic cycle L a U L b U T ab U T ba corresponds to a 
uniquely defined sequence of points Mj £ S a U £& such that 

M i+1 = 

where 

6 = c if Mj e S c (ce{a,b}). 

In this way the trajectory of the frozen system defines a sequence {&K=-^ which 
is called the code of the orbit. 

The periodic orbits L a and L b are saddle and the intersections of the stable and 
unstable manifolds of L a and L b that create the heteroclinic orbits are transverse 
due to the assumptions [Al] and [A2]. Consequently (cf. [3j), one can choose the 
maps Yi ab and U ba and define coordinates (x, y, z) in S a and T, b in such a way 
that the following holds. 

[HI] S c fl M.q is diffeomorphic to the product X c x Y c x D where X c and Y c are 
balls in M m_1 of a radius R > centred around the origin. 

[H2] For each pair c, c' £ { a, b } the Poincare map U cc r can be written in the 
following "cross-form" [12] : there exist smooth functions f cc i , g cc i : X c x 
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Yd — > X c / x F c such that a point M(x, y, z) G S c is mapped to M(u, w, z) G 
E c / by the map II CC / if and only if 

x = fcc'(x,y,z), y = 9ccf(x,y,z). (18) 
[H3] There exists A < 1 such that 



d(fce> 


9cd) 


d(x 


y) 



where the norm of the Jacobian matrix corresponds to max{||x||, ||y||}. 

Inequality f|T9|) implies that for a fixed z G D the set of all orbits that lie 
entirely in a sufficiently small neighbourhood of the heteroclinic cycle L a U Lb U 
r a b U Tba in the energy level A4q is hyperbolic, a horseshoe. Moreover, one can 
show that the orbits in A z are in one-to-one correspondence with the set of all 
sequences of a's and b's, i.e. for every sequence {£1}]=!^ there exists one and 
only one orbit in A z which has this sequence as its code. 

Indeed, take any orbit from A z and denote by Mj(xj, y^ z) the sequence of its 
intersections with the cross sections. Equation f)18p . implies that the orbit has a 
code {CiYiZ-x, ^ an d only if the coordinates of satisfy the equations 

x i+i = ftei+i( x i,yi+iiz), v% = 9aA +1 ( x iiyi+ii z ) ■ 

Therefore the sequence { (xj,yj) j^oo is a fixed point of the operator 

{(zt>3/i)}£=-oo i - > {(hi-^i( x i-uy^ z ),9^ i+1 { x i,yi+i^)}tZoo- ( 20 ) 

Equation (|T9|) implies this operator is a contraction of the space n J^oo x ^& > 
hence the existence and uniqueness of the orbit with the code {£1}]=!^ follow 
(see e.g. [IlJ). 

Moreover, the fixed point of a smooth contracting map depends smoothly on 
parameters. Consequently the orbit depends smoothly on z and the derivatives 
of (xi(z,£),yi(z,£)) are bounded uniformly for all i. 

Lemma 2 If the Poincare maps satisfy assumptions [H1]-[H3], then for any two 
code sequences ^ = }t^-oo an< ^ £^ = {£i^}i=^oo> which satisfy 

Q Q J or M — n 

the corresponding intersections with the cross sections are bounded by 

max{ ||*^,e (1) ) -Xi(z,^)\\, \\yi(z,^)- yi (z,^)\\ } < 2RX n ^ , (21) 

where the constants R > and A G (0, 1) are defined in [HI] and [H3] respectively 
and do not depend on the sequences f;^' 2 '. 
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Proof. First we note, that 

max { \\xi(z, £ {1) ) - x t (z, £<*>) || , \\y i+1 (z, - y i+1 (z, } 
< Amax{||x m (2,e (1) )-x m (z,e (2) )||, \\y t (z, £«) - y t (z, £ (2) ) || } (22) 

for — n < i < n. Since none of the normes involved exceeds 2R we immediately 
conclude that 

max{ \\x t (z,^) - x 4 (z,e (2) )||, \\y i+1 (z, £ (1) ) - y i+1 (z, £ (2) ) || } < 2/?A . (23) 
Then the following estimate is true for n' = 1 

max{ -Xi(z,t&>)\\, \\y i+ i{z,ZW) - y i+l {z^ 2) )\\ } 

f A n ' _ * for < ?' < r?' 
<2i? , "/ i • (24 ) 

[ A n for —n' - 1 < 2 < v 7 

We continue inductively in n'. Assuming that the estimate (1241) holds for n! 
replaced by n' — 1 we check the upper bounds for all 2n' + 1 different values of 
i in the order of decreasing of On each step we use the contraction property 
( 122]) and the sharper of upper bounds ( 124|) and ( 1231) . In the case n 1 = n, 

max{ \\x n (z,e ] ) - x n (z,^)\\, h-niz,^) - V-n(z, e 2) ) II } < 2R 



is used instead of ( 1231 . Then (l2lj) follows directly from the last upper bound and 
(1241) taken with n' = n. □ 



Now let us consider the full system (CQ) for a small e > 0. Since the vector 
field depends smoothly on e, the Poincare maps U cc r : S c — > S c / are still defined 
and can be written in the following form: 

x = fcc'(x,y,z,e), y = g cc >(x,y,z,e) 
z = z + E(p cc >{x,y,z,e), 

where f,g,<fi are bounded along with the first derivatives and f,g satisfy (1TU1) . 

For technical reasons we need to assume that the domain D is invariant under 
the Poincare map, i.e., cc '(x, y, z, e) = if z G 9£). If this is not the case, then 
we modify <p cc > in a small neighbourhood of the boundary. We note that the 
next lemma contains a statement of uniqueness but the surfaces provided by the 
lemma may depend on the way the functions (p cc i have been modified. Therefore 
the lemma implies existence but not uniqueness for the original system. 

The next lemma is of a general nature and has little to do with the Hamilto- 
nian structure of the equations. Rather we notice that by fixing any code £ and 
varying z G D we obtain at e = a sequence of smooth two-dimensional surfaces. 
The i-th surface is the set run by the point Mj(z) of the uniquely defined orbit 
with the code £. This sequence is invariant with respect to the corresponding 
Poincare maps and is uniformly normally-hyperbolic — hence it persists at all e 
sufficiently small. 
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Lemma 3 Given any sequence £ of a's and b's, there exists a uniquely defined 
sequence of smooth surfaces 

e) : [x, y) = (xi(z, £, s),yi(z, £, e)) (26) 

such that 

n^Ci+iA = A+i- (27) 

T/ie functions (xi, yi) are defined for all small e and all z G D, they are uniformly 
bounded along with their derivatives with respect to z and satisfy (I2~TT) . Moreover, 
there is C > independent from the code £ such that 

\\(xi(z,€,e) -Xi(z,£,0) t yi(z,€,e) - y l {z,^,0))\\ < Ce , 

for all i G Z. 

A proof of this lemma is essentially identical to the proof of Lemma 1 of [6J and 
is based on contraction mapping arguments: the functions Xi,yi are constructed 
as a fixed point of an operator similar to (1201) . 

We note that if £ = c°° is a code which consists of the symbol c only, then 
Ci, Xi and y^ are independent from i, and we will denote them by £ c , x c and y c 
respectively 

5 Drift of slow variables 

Let £ be a code. The corresponding trajectory of the full system is described by 
the dynamics of its slow component: 

Zi+i = Zi + e(j} iiZi+x {x i {z i ^,e),y i {zi^,e),z h e) . (28) 

If £ = c°° the functions x% and yi do not depend on i, so we denoted them by x c 
and y c respectively. Then the equation can be written in the form 

z i+ i = Zi + e(f) cc (x c (z u e),y c (zi, e), z h e) , (29) 

where the bars over Zi and z i+1 are used to distinguish trajectories of (129]) and 

(EHD- 

The next lemma estimates the difference between these two slow dynamics 
for all sequences which have a large block of c's. 

Lemma 4 Assume the assumptions of Lemma are satisfied. Then for any 
K Q > 0, t > 0, there is K\ > and e > such that for any \e\ < eo and any 
code £ such that for some index j 

= 0+i = • • • = £j+L^j = c 
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the inequality \\zj — Zq\\ < sKq implies the corresponding trajectories of (|28|) and 
(|29|) satisfy the inequality 

\\zj+N ~ zn\\ < eKi 

for allO<N< N (e) = [|J . 
Proof. Using (|28|) we get 

j+N-l 

Z j +N = Z j +e ^ <l>tei + i( x i( z ht> £ ),yi( z ii€, £ )i z U £ ) ■ (3°) 



Using we obtain in a similar way 



N-l 



(31) 



i=0 



We have assumed \\z — Zj-|| < K e. Then taking the difference of the equalities 
f l30|) and fT3Tj) . we obtain 



iV-l 



\Zj+N ~ Z N\ 



< eK + e ^2\\<f> cc{,-Ej+i(.Zj+ii £J j Vj+i^Zj+i, s), ^j'-fi; £J 



Consequently, 

\\z j+N -z N \\ < eK + e 

d<j> a 



dz 



i=0 

-0 cc (x c (^, e), yJzi, e),Zi, e) 



N-l 

y^, \\zj+i - Zi 



(32) 



(33) 



i=0 



+e 



d{x,y) 



N-l 



II ( x j+i( z j+i> £ ) — x c(zi, e), yj + i(zj +i , £, e) — y c (zi, e) 



i=0 



In order to estimate the last term we note that Lemma[3]includes the estimate (I2T1) 



JV-l 



8=0 



^] II (^i+K^i+i) £ ) ~ ^(-^ii £ )i Uj+i(zj+i, £, e) — y c ( z i, e))\\ 
=o 

— II £> £ ) — ^c(^j+i) £ )i Uj+i(Zj+i, £, e) — y c (zj + i, e) 



i=0 
N-l 



+ ^2\\( x c(z j+i ,e) 



x c {Zi, £), y c {Zj+i, £) — y c {Zi, £ 

i=0 

N-l N-l 

< 2RX min{i ' Noi£) - i} + max 

i=0 i=0 



/ 4i? 

< + max 

1 — A 



dz 



dy c 



dz 



N-l 

E 

i=0 



Zj-~i | 
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Substituting the last bound into (1331) we obtain 
||%+iV — Zn II < £ 





d(j) C c 


+ 




max | 


dx c 




dy c 


} 


( 


dz 




d(x,y) 


dz 




dz 





N-l 

£ 



I 



+£ 



a(x,?/) 



AR 
1- A 



Let 



.4 
5 



<9(x,y) 



4f? 
1- A 



<9</> cc 




d(j) cc 


max | 


<9x c 






) 


dz 


+ 


d{x,y) 




• 


(9,2 





Then 



JV-l 



\z j+N - z N \\ < eA + ^ 



i=0 



and consequently 

So for iV < ^ we have 

— e 

where A4 = Ae*° B . 



\zj +N ~ z N \\ < eAe 
\\zj+N ~ z N \\ < K x e 



(34) 
□ 



We note that Lemma H] is also valid for any two sequences with any large 
common block. 

Lemma 5 Assume the assumptions of Lemma are satisfied. Then for any 
Kq > 0, to > 0, there is K\ > and Eq > such that for any \e\ < £q and any 
two codes ^ and ^ such that for some index j 



ejX = ejl o < i < n (e) 



to 

E 



the inequality 



7 « _ r ( 2 ) 



< eK implies 



for allO<N< N {e). 

The proof of this lemma is almost identical to the previous one so we skip it. 
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6 Proof of Theorem H 



Now we have all ingredients necessary to complete the proof of Theorem [TJ Each 
periodic orbit L k E A is defined by a periodic code. Let £o be the longest period 
among the codes corresponding to the periodic orbits selected in the assumptions 
of Theorem [TJ Then each of the periodic orbits can be uniquely identified by a 
piece of code Ck E {a,b } e °. 

Given an accessible path F we define 

Aj = Tj+i - Ti . 

It is the time the slow motion follows the flow defined by the Hamiltonian function 
Jk t i 1 < i < N, where A" is the number of segments in the path. Let 



Ni{e) 



A,; 



el 



Let u>i = be a finite sequence, which consists of iV»(e) copies of the symbol 
Cfc.. Let £ £ be any sequence, which contains starting from position 0. 

Obviously, the sequence 

i 

1=0 

indicates the starting positions of the blocks u>i in the code £ £ . 

We note that assumptions [Al] and [A2] imply that there are sections and 
Poincare maps of the frozen system (J2]) which satisfy [HI] — [H3]. In order to apply 
Lemma [3] we have to modify the slow component of the Poincare maps to ensure 
invariance of D. Since D is open there is 5 > such that a ^-neighbourhood of 
T is inside D. Then we modify <p cc ' outside this 5-neighbourhood of F to ensure 
that <p cc i vanishes near dD. This modification does not affect trajectories which 
do not leave an 0(e) neighbourhood of F: i.e. while a trajectory of the modified 
maps stays in the neighbourhood of F it is simultaneously a trajectory of the 
original Poincare maps. 

Lemma [3] implies that there is a sequence of surfaces which corresponds to 
the sequence £ = £ e . Now let z = T(0) and consider the sequence of points 

Mi = M^{xi{zi,i,e),yi{zi,^e),Zi\e) 

on those surfaces. The slow component Z{ satisfies equation ( 1281) and 

Mi E M £ fl E 6 . 

Lemma H] (or Lemma [5]) implies that 

\\zi - Zi\\ < K x e 
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for all i such that < % < Ni(e), where Z{ denote the trajectory of (1291) . We 
continue inductively to show using Lemma H] that there are constant such 
that 

IN - ^-ifc-ill < K k£ (35) 
for all % such that jk-i < i < jfc and 1 < A; < A/", where Zi satisfy (1281) with 
initial condition := ^i fe _i- We note that these Zi-j k _ 1 all belong to the invariant 
surface £ Cfc , then Lemma [3] implies 

dist(M i ,5 Cfci0 ) = 0(e). 

We consider the trajectory of the full system ([I]), which we denote by j £ , such 
that 7 e (0) = M . Since 7 e goes through the points Mj it also stays 0(e)-close to 
S Ck fi between the points therefore 

diBt(7 e (*),5 Ch)0 ) = O(e) 

for t G [Tk-i£~ l , TkS' 1 }. Then Lemma [T] implies that the slow component n s 'y £ (t) 
shadows the accessible path T. 



7 Final remarks 

Finally, we note that similar equations arise in the case of a Hamiltonian system 
with the standard symplectic form, Q st = dp A dq + dv A du when a Hamilto- 
nian function looses some degrees of freedom at e = 0. More precisely, if the 
Hamiltonian function has the form 

H(p,q,v,u; e) = H (p,q) + £H x {p, q, v, u) , 

the corresponding Hamilton equations are given by 

dH dHx dH dH x 

dp dp dq dq ,o fi \ 

dH x dH, W 

dv du 

In this equation, adiabatic invariants can be destructed by resonances [10] . 

These equation are quite similar to ([T|). We note that the frozen fast system 
is independent of the slow variables. The theory developed in this paper can 
be applied to the system (1561) . The most notable difference is related to the 
description of the slow motion near a cylinder formed by periodic orbits of the 
frozen system. Indeed the slow motion is described by the averaged perturbation 
term 

J c (u,v) = Hi{p c (t),q c (t),u,v)dt 
Jo 

and not by the actions. 
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